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Abstract. Quantum algorithms for several problems in graph theory are consid- 
ered. Classical algorithms for finding the lowest weight path between two points in 
a graph and for finding a minimal weight spanning tree involve searching over some 
space. Modification of classical algorithms due to Dijkstra and Prim allows quantum 
search to replace classical search and leads to more efficient algorithms. In the case of 
highly asymmetric complete bipartite graphs, simply replacing classical search with 
quantum search leads to a faster quantum algorithm. A fast quantum algorithm for 
computing the diameter of a complete graph is also given. 



Introduction 

The question of which classical algorithms can be sped up by quantum com- 
puting is of course a very interesting one. At present there are only a few general 
techniques known in the field of quantum computing and finding new problems 
that are amenable to quantum speedups is a high priority. Classically, one area of 
mathematics that is full of interesting algorithms is computational graph theory. 
It is therefore natural to ask whether any of the classical graph theory algorithms 
can take advantage of quantum computing. 

One of the few general techniques known centers around Grover's algorithm for 
searching an unsorted list for a specified element. This original idea has been 
extended to general amplitude amplification that can be applied to any classical 
algorithm. It would be incorrect to assume that amplitude amplification always 
leads to quantum speedups of classical algorithms. There are some interesting cases 
where "Grover-like" techniques do in fact lead to speedups of classical algorithms. 
One very important case of this is to find the minimum value of a computable 
function as the set of input arguments ranges over a finite, but unordered list. In 
this case, if the list is of length n, then the quantum cost of finding the minimum is 
0[\/n), while the classical cost is 0(n). Quantum algorithms for searching for the 
maximum or minimum of an unsorted list have been described in [DH] and [AK]. 
The question to be addressed here is whether this leads to speedups in graph theory 
algorithms that employ classical minimum finding in the course of solving a graph 
theory problem. 
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1. Minimal Weight Paths 

The problem of finding the shortest path between two points in a weighted 
graph is an old one. If G = {V, E) is a directed graph with a weighting function 
v : E ]R+, the weight of a path is the sum of the weights of the edges that 
comprise the path. If G is a complete graph, then the function i/ is well defined on 
all the edges. In case G is not a complete graph, and {v,v') ^ E, then it is useful 
to define v{v, v') = oo. This allows u to be defined on all ofVxV, not just E. 

Pick 3. vq E V from which all the shortest paths are to be computed. Consider 
the following algorithm due to Dijkstra (see [G]): 

(1) S ^ {vo}, X{vo) ^ 

(2) for u e F - 5 do 
(2.a) X{v) ^ iy{vo,v) 

(3) while S do 

(3. a) find w eV — S such that X{w) is minimal 
(3.h) S^Syj{w} 
(3.c) for we y-S' do 

(3.C.1) \{v) ^ min(A(i;), A(m;) + y{w,v)) 

Figure 1: Dijkstra' s algorithm. 

At the end of this procedure \{v) is the length of the shortest path from to v, and 
with only minor changes, this algorithm is easily modified to record the information 
needed to construct the shortest path. This algorithm consists of iterating over all 
of V by successively adding elements to S, which is the set of points for which 
the shortest path from fo has already been determined. Each iteration consists 
of a search procedure to find the next nearest element of to wq, and an update 
procedure to record all the newest shortest path information for the remaining 
vertices based on this newest next nearest clement. 

Analysis of the run time for this algorithm is quite simple. At line (3) if z = 15 |, 
then the (classical) cost of the search in line (3a) is n — i and the cost of updating 
in line (3c) is n — i — 1. The total cost is therefore Y11=q 2n — 2z — 1 = O(n^) . 

There are some important modifications to this original algorithm of Dijkstra if 
the graph being searched is somewhat sparse (i.e. if the total number of edges is 
much less than n^). In this case, Dijkstra's algorithm can be modified to find the 
shortest path with work 0((|y | + |i? |) log \V \ ) through the use of priority queues. 
In general, the key idea is that the update procedure in line (3c) need only update 
the shortest path for those vertices that arc adjoining the vertex w which was most 
recently added to the set S of those vertices whose shortest path distance from vq 
have alread been computed. The reason the only the vertices adjacent to w need 
to be considered is that for all other vertices j/(w,f) = oo. 

Unfortunately, if the classical search for the minimum at line (3a) is replaced by 
the quantum algorithm for finding the minimum of an unordered set, the cost for 
line (3a) per iteration drops to 0[\/n — , but the cost of the entire algorithm is 
still O(n^) since the update cost per iteration in line (3c) is still 0{n — i). 

One possible way around this problem is to dispense entirely with the update 
procedure in line (3c) at the cost of a larger search in line (3a). This modified 
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algorithm then goes as follows: 



(1) S ^ {vo}, X{vo) ^ 

(2) while S^V do 

(2. a) find {w,v) G S* x {V — S) such that \{w) + v{w,v) is minimal 
(2.h) 5 5 U {v}, \{v) \{w) + v{w, v) 



As before, at the end of this procedure \{v) is the length of the shortest path 
from vq to and again, this algorithm is easily modified to record the information 
needed to construct the shortest path. 

The analysis of this modified algorithm is also quite easy. At line (2) if i = \S\, 
then the (classical) cost of the search in line (2a) is i {n — i), making the entire 
classical cost of this algorithm X]r=o iin — i) = 0(n^), which is quite a bit worse 
than Dijkstra's original algorithm. However, the quantum cost of this algorithm is 
determined by noting that the search cost at line (2a) reduced to ©(i/i (n — i)), 
thereby making the total cost of the algorithm O (^"^g — i)) = O(n^) , which 

can be seen by noting that 



doesn't represent an improvement over the original classical algorithm. 

What seems to be really needed is to have a partial tradcoif between the search 
and update parts of the algorithm. The idea is to balance the classical update cost 
per iteration with the quantum search cost. The following algorithm is one way of 
accomplishing this. 

(1) S ^ {vo}, T ^ {vo}, X{vo) ^ 

(2) for w e V - 5 do 
(2.a) X{v) ^ iy{vo,v) 

(3) while S^V do 

(3. a) find (?i7,w) G T x [V — S) such that \{w) + v{w,v) is minimal 
(3.b) find u & V — S such that A(u) is minimal 
(3.c) if \{w) + i'{w,v) < A(w) then 

(3.C.1) S ^ SU{v}, T ^TLl{v}, X{v) ^ X{w) + iy{w,v) 
(3.d) if X{w) + v{w,v) > X{u) then 
(S.d.l) S^SU{u}, T^TU{u} 
(3.e) if |T| > A; then do 
(3.e.l) for w e V - S' do 

(3.e.l.a) find w G T such that X{w) + v{w,v) is minimal 
(3.e.l.h) if X{w)+v{w,v) < X(v) then 
(3.e.l.b.l) X{v) ^ X(w) + v{w,v) 
(3.e.2) T ^ {vo} 

Figure 3: Dijkstra's algorithm with periodic updating 



Figure 2: Dijkstra's algorithm without full updating 
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The idea is to have a set T of vertices for which the full update of X{v) for all 
remaining vertices inV — S has not yet been computed. This full updating is done 
every k-th iteration of the main loop. By keeping vo in T all the time, there is 
always the possibility of going directly from vq to v since X{vo) = 0. The value of 
fc is a parameter for this algorithm and needs to be set to optimize the total cost. 

Most of the work in this algorithm takes place in line (3), the main iteration, 
which is done a total of n times. To analyze the work for the i-th iteration, write 
i = hk + j with 1 < j < A;, so that the size of the set T on the i-th iteration is j. 
It is convenient for this analysis to assume that k divides n since that makes line 
(3c) execute exactly n/k times, however, even if this assumption does not hold, the 
work calculation is still valid. 

The search for the minimum on line (3a) is over a set of size j {n — i), so the 
total work over all iteration of line (3a) is 

n/fc-l fc 

^ J2j{n-hk-j)=0{kn^) 

h=0 j=l 

in the classical case and 

n/fe-l k 

E T.Vj{n-hk-j) = 0{k'/'n^/') 

h=0 j=l 

in the quantum case. 

The update cost on line (3el) requires n — i searches for the minimum over a 

set of size k each time and this is done only when i is divisible by k. Therefore 
the entire cost of the update procedure on line (3e) over all the iterations of the 
algorithm is 

n/k 

^{n-kh)k = 0{n'^) 

h=l 

in the classical case and 

n/k 

j2{n-kh)Vk = 0{k-^^^n'') 

in the quantum case. 

In the classical case, the total work for the algorithm is max(0(fcn^), 0(n^)) = 
O(n^) which is minimized by taking fc = 1, and therefore 

Wclassical = O(n^). 

Note that taking fc = 1 in figure 3 gives the original Dijkstra algorithm of figure 
1, while taking fc = |F| yields the algorithm of figure 2, so fc might reasonably be 
viewed as an interpolation parameter. 

In the quantum case, the situation is a bit different. The total work in line (3) 
is just the maximum of the work in lines (3a) and (3e), since the work in line (3b) 
is always dominated by these other work factors. The total work is therefore 

max(0(fci/2n-''/2),0(fc-i/2„2)) 

and to minimize this, the parameter fc should be chosen to make these two work 
factors the same. Setting fc^/^n^/^ = fc~^/^n^ gives fc = n^/^ and therefore 

This indeed is an improvement over the classical work factor of n^. 
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2. Graph Diameter 

The algorithms so far presented compute the minimal path length from one point 
in the graph to another point in the graph. This minimal path length will be referred 
to as the distance from one point to another in the graph. The diameter of a graph 
is the distance bc;twc;c!n the two furthest points in a graph (i.e. the maximum of the 
minimal path lengths in the graph). This section deals with quantum algorithms 
for finding the diameter of a weighted graph. 

Since the initial point vq was always fixed in Dijkstra's algorithm and its variants, 
this was never explicitly indicated in the notation X{v). However for the purposed of 
this section, it is useful to write X{vo, v) in place of X{v). One of the key features of 
all of the "Dijkstra-like" algorithms described in the first section is that they start 
from a given vertex vq and by an iterative procedure manage to find the minimal 
distances to all the other vertices in the graph. Furthermore the last vertex for 
which the distance is computed is always the most distant vertex from vq. These 
"Dijkstra-like" algorithms could therefore be viewed as computing the maximum 
distance from vq- The diameter of the graph is just the maximum of all these 
distances as vq runs over all of V. By invoking the quantum maximum finding 
algorithm with a "Dijkstra-like" algorithm as a callable subroutine, the diameter 
of the graph will follow. Since there are n = |F | possible initial values of vq, the 
quantum cost is simply ^/n, times the cost of the inner loop. Using the best quantum 
algorithm as the inner loop gives a total quantum cost for finding the diameter as 
0(n9/4). 

By way of comparison with classical costs, Dijkstra's (classical) algorithm would 
have to be run n times giving a classical cost of O(n^). There is an interesting 
alternative classical algorithm due to Floyd and Warshall that finds all the distances 
between all pairs of vertices in a weighted graph. Its cost is also O(n^). 

3. Minimal Weight Spanning Trees 

Another common problem in classical graph theory is that of finding a minimal 

weight spanning tree of a graph. There is a very nice classical algorithm for this 
due to Prim that goes as follows (see [G]): 

(1) S ^ {v,}, F^{} 

(2) iovv^V -S do 

(2.a) L{v) ^ iy{vo, v), M{v) ^ vq 

(3) while St^V do 

(3. a) find w gV — S such that L{w) is minimal 
(3.b) S^SLl{w}, F ^ Fu{{w,M{w)} 
(3.c) for G V - S' do 

(3.C.1) if p{w,v) < L{v) then 

(3.c.l.a) L{v) i>{w,v), M{v) <- w 

Figure J^: Prim's algorithm. 

At the end of this algorithm, f is a set of edges that comprise a minimal weight 
spanning tree of the original graph {V, E), which has been implicitly assumed to be 
a complete graph on n = \ V\ vertices. As in the case of Dijkstra's minimal weight 
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path algorithm, this algorithm consists of an scarc;h phase for each iteration and 
an update phase for each iteration. In the original algorithm of Prim, just as in 
Dijkstra's original algorithm, the work for the search phase and the work for the 
update phase are carefully balanced to be the same, the total work being O(n^). 

There is a version of this algorithm that doesn't employ any updating at the 
cost of having to search a larger space to find a minimum on each iteration. That 
algorithm goes as follows: 

(1) S ^ {v^}, F^{} 

(2) while 5 7^ V do 

(2. a) find {u,v) £ S x {V — S) such that i'{u,v) is minimal 
(2.b) S^SU {v}, F^FU {{u, v)} 

Figure 5: Prim's algorithm with no updating 

Essentially the same idea that works for a good quantum version of Dijkstra's 
algorithm, as in the first section of this note, works to give a good quantum version 
of Prim's algorithm. In particular, the following modification of Prim's algorithm 
is what seems to work: 

(1) S ^ {vo}, T ^ {vo}, F ^ {}, L{vo) ^ 

(2) ior veV-S do 

(2. a) \{v) ^ v{vo,v), M{v) ^ vq 

(3) while 5 ^ F do 

(3. a) find G T x [V — S) such that v{u,v) is minimal 

(3.h) find w €V — S such that L(w) is minimal 
(3.c) \f L{w) < v{u,v) then 

(3.C.1) Syj{w], T<-TU{w}, F ^ Fyj{{w,M{w)} 
(3.d) if L{w) > z/(u, w) then 

(S.d.l) S\J{v}, T^T\j{v}, F^F\j{{u,v)} 
(3.e) if |T| > fc then do 
(3.e.l) for V eV S do 

(3.e.l.a) find w such that i'{w,v) is minimal 
(3.e.l.h) if v{w,v) < L{y) then 

(3.e.l.h.l) L{v) <- v{w,v), M{v) w 
(3.e.2) T ^ {vq} 

Figure 6: Prim 's algorithm with periodic updating 

The analysis of this algorithm is practically identical to that of Dijkstra's algorithm, 
and again the parameter k is determined to equalize the work of searching for a 
minimum with the work of updating. If quantum minimum finding is used in steps 
(3a), (3b), and (3ela) instead of classical minimum finding, the optimal value of k 
is again 0(n^/^) and the overall work for the quantum algorithm in figure 6 here 
is 0(n^/*), which is an improvement on the classical algorithm. 
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4. Bipartite Graphs 

Another approach is to change the problem to one that is more amenable to 
replacing the classical search with a quantum search. The idea is to balance the 
classical update cost per iteration with the quantum search cost by introducing 
asymmetry into the problem in Dijkstra's original algorithm. 

In a bipartite graph, the set of vertices V is divided into two disjoint subsets Vi 
and V2, and all edges connect points in Vi with points in V2 (and vice versa, as well 
in a directed graph). In a complete bipartite graph, each point in Vi is connected 
to every point in V2, and conversely, each point in V2 is connected to every point 
in Fi. If rij = |Vi I for i = 1, 2, then the number of edges in a complete bipartite 
graph is nin2, and in a complete bipartite digraph, the number of edges is 2nin2- 
If rti and n2 are of very different sizes, then a complete bipartite graph will be quite 
sparse. In what follows, for this section, the question of finding the lowest weight 
path in complete bipartite graphs will be considered. Since a complete bipartite 
graph has edge set = (Vi x V2) U (V2 x Vi), it will be assumed that a pair of 
weight functions ui : Vi x V2 ^ K"*" and U2 : V2 x Vi ^ R+ are given. The length 
of the shortest path Xi : Vi ^ R+ from an initial vertex vq to elements in Vi for 
i = 1,2 will be computed by the algorithm. 

Assuming that ni < n2, the work for Dijkstra's algorithm is 0(nin2 logn2), 
with a sparse graph modification of the algorithm using priority queues. There 
are certain cases in which this may not be the best classical algorithm for finding 
the shortest path between two points. The run time for Dijkstra's algorithm may 
be optimal, but the memory requirements may be overwhelming. The following 
algorithm constructs the list of shortest distances (and paths if a little additional 
information is kept) between a fixed initial point vq S Vi and all the other points 
in Vi U V2- 

(1) S^^b^'o}, Xi{vo)^0 

(2) for w e Vi - 5i do 

(2. a) find w e V2 such that ui{vq,w) + V2{w,v) is minimal 
(2.h) \i{v) ^ vi{va,w) + V2{w,v) 

(3) v/hi\e Si ^ Vi do 

(3. a) find {u, v, w) e x V2 x {Vi-Si) such that \i{u) + ui{u,v) + V2{v,'w) 
is minimal 

(3.b) 5i ^ 5i U {w}, \i{w) ^ \i{u) + vi{u,v) + V2{v,w) 

(3.c) for t; e Vi - 5 do 

(3.C.1) find u e V2 such that Xi{w) + ui{w,u) + U2{u,v) is minimal 
(3.C.2) Ai(w) ^ vam{\i{v),Xi{w) + vi{w,u) +V2{u,v)) 
U) for ue ^2 do 

(i-a) find w & V\ such that Ai(«;) + vi{w, u) is minimal 

(4.h) \2{u) <— Xi{w) + ui{w,u) 

Figure 7: Dijkstra 's algorithm for a bipartite graph with partial updating. 

The idea is that since Vi is smaller than V2, only the distances in Vi need to be 
updated. Vertices in V2 are viewed merely as intermediate points along the way for 
paths that connect vertices in Vi. 
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The analysis of this algorithm is only slightly more complicated than the previous 
analyses. Beginning with the idea that classical search for the minimum of some 
function over a set of size N requires 0[N) operations, it is clear that line (2a) 
requires work 0(712) for each iteration of line (2), and since there are ni iterations of 
line (2), the total work for the initialization of X{v) for all v &Vi, which is what goes 
on in lines (1), (2), (2a), and (2b), is 0(nin2). As for line (3), let i = \S \ on each 
iteration. Then for (3a), the search is over a set of size i (ni — i)n2, which makes 
the total search cost for (3a) over all the iterations in line (3) '^^Iq i (ni — i)n2 = 
0(nfn2). As for the update operations that start in line (3), there is an subitcration 
over a set of size ni — i and within each subiteration, there is a search over a set 
of size n2- Thus the total update cost is also J27=o ^ ("1 ^ "^2 = 0{n\n2) ■ Thus 
the total (classical) algorithm cost is 0{n\n2), but only a memory of size 0(ni) 
is needed. For the final stage (4) of the algorithm that fills in the cost function 
for the elements of the second and larger part, V2, the outer iteration is over a set 
of size 712. For each u G V2, the minimum in line (4a) is found over a set of size 
ni, making the total (classical) work for this final phase of the algorithm O (711712), 
which is clearly not the dominant cost. 

Now for the quantum costs. The costs of initialization are ni iterations of search 
over a set of size 712, so the quantum cost is 0{n\n^^^) . The search cost per iteration 
in line (3a) is 0{^\/i {n\ — i) 712), making the total search cost 

V i {ni — i) 7i2 » 7i2^^ / ^yx{n — x)dx = O^nfn^^^) . 

The update procedure, starting on line (3c) consists of subiterating over a set of 
size ni — i, and for each subiteration, there is a search over a set of size 712- The 
total cost of updating is therefore 

ni i 

i=o i=o 

which is the same as the search cost. The total cost of this algorithm prior to 
computing the minimum weight paths for all of V2 is therefore 0(n\ny^), with 
a memory cost of 0(7ii). The final stage of the algorithm costs 712 iterations of 
minimum finding over a set of size 711, so the (quantum) work for filling in the cost 

function for V2 is 0(71271^^). The total cost of the whole algorithm is therefore 

r\i 2 1/2 , l/2\ 
U{n{n2 +71271]^ j. 

How does this compare to the best classical cost, which is 0(711712) with a mem- 
ory cost of 0(711772), as well? Suppose that only the minimum weight paths to the 
other elements of Vi are desired. Then if 0(77? 77^^) < 0( 771772), then the quantum 
algorithm wins. This occurs if Tii = o(y^). Thus for highly unbalanced complete 
bipartite graphs, the quantum algorithm outperforms the classical algorithm. 

Conclusion 

The quantum versions of Dijkstra's algorithm in section 1 and the quantum ver- 
sion of Prim's algorithms in section 3 are for complete graphs with a well defined 
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weight function on all the edges. The improvements over the corresponding classi- 
cal algorithms are also for complete graphs. If a graph is not complete, then the 
classical minimal weight path algorithm and the classical minimal weight spanning 
tree algorithms can effectively use the graph structure to give better classical algo- 
rithms. Essentially this is done through updating only the vertices that are adjacent 
to the vertex being added to the core set. Classically, the use of priority queues lead 
to these improvements. It seems likely that these ideas can be incorporated into 
future quantum algorithms for incomplete graphs, but this is a topic for a future 
paper. One particular case of this for complete bipartite graphs was dealt with in 
section 4. 

There are a number of closely related problems to finding minimal weight paths 
in graphs. Among these are the problems of deciding graph connectivity and finding 
the shortest path between two points. Determining the number of blocks in a graph 
follows directly from being able to decide graph connectivity as does determining 
the set of articulation points in a graph. Presumably the quantum algorithms given 
here can be readily extended to these problems, although the details are yet to be 
worked out. 

The whole field of quantum algorithmic graph theory has barely been touched on 
here. However, the algorithms developed in this note should be regarded as strong 
evidence that many well known classical graph theory algorithms have interesting 
quantum analogues. Often they will involve some nontrivial modifications of the 
classical algorithm to make optimal use of the few tools currently available in the 
quantum toolbox. The optimistic view is that although the modifications to the 
classical algorithms may be nontrivial, they often are not excessively complicated 
either, as in the case of the quantum versions of Dijkstra's algorithm and Prim's 
algorithm. 
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